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Abstract 

Quantum-correlated p(3770) — DD decays collected by the GLEO-c experi¬ 
ment are used to perform a first measurement of the fractional CP-even 
content of the self-conjugate decay D —)■ 7r+7r“7r+7r“, obtaining a value of 
0.737 ± 0.028. An important input to the measurement comes from the use 
of P —)■ A'g 7 r’'' 7 r“ and D —)■ K^n^n~ decays to tag the signal mode. This 
same technique is applied to the channels D —)■ 7r’''7r“7r° and D —)■ 
yielding P^^^° = 1.014 ± 0.045 ± 0.022 and = 0.734 ± 0.106 ± 0.054, 

where the first uncertainty is statistical and the second systematic. These 
measurements are consistent with those of an earlier analysis, based on CP- 
eigenstate tags, and can be combined to give values of = 0.973 ± 0.017 
and = 0.732 ± 0.055. The results will enable the three modes to be 

included in a model-independent manner in measurements of the unitarity 
triangle angle 7 using —)■ DK^ decays, and in time-dependent studies of 

CP violation and mixing in the system. 
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1. Introduction 


Studies of the process —)■ DK^, where D indicates a neutral charmed 

meson reconstructed in a state accessible to both and decays, give 
sensitivity to the unitarity triangle angle 7 = arg(— 14 idKib/^dKb) (also de¬ 
noted Improved knowledge of 7 is necessary for testing the Standard 
Model description of CP violation. In a recent publication [1] it was shown 
how inclusive three-body self-conjugate D meson decays can be used for this 
purpose, provided their fractional CP-even content is known, a quantity de¬ 
noted P+ (or fI when it is necessary to designate the specihc decay /). 
Measurements of P+ for the decays D —)■ 7 r+ 7 r“ 7 r‘’ and D —)■ were 

performed, making use of quantum-correlated DD decays coherently pro¬ 
duced at the 'ip (3770) resonance and collected by the CLEO-c detector. In 
this Letter a hrst measurement is presented of the CP content of the four- 
body mode D —)■ 7 r+ 7 r“ 7 r+ 7 r“, again exploiting CLEO-c ^’(3770) data. This 
fully-charged and relatively abundant hnal state [ 2 ] can be reconstructed 
with good efficiency by the LHCb detector and hence is a promising mode 
for improving the determinaton of 7 at that experiment, as well as at Belle 11. 

The three-body analysis reported in Ref. |1] exploited events in which 
one D meson is reconstructed in the signal mode and the other ‘tagging’ me¬ 
son in its decay to a CP eigenstate. The measurement of presented 
in this Letter follows the same method, but augments it with other ap¬ 
proaches, in particular a complementary strategy in which the tagging modes 
are D —)• and attention is paid to where on the Dalitz plot this tag 

decay occurs. In order to beneht from this strategy for the previously stud¬ 
ied decays, this Letter also presents measurements of and using 

D Kl tags. Throughout the effects of CP violation in the charm 

system are neglected, which is a good assumption given theoretical expecta¬ 
tions and current experimental limits Pjl] . However, as discussed in Ref. |5] , 
knowledge of P+ also allows such D decays to be used to study CP-violating 
observables and mixing parameters through time-dependent measurements 
at facilities where the mesons are produced incoherently. 

The remainder of the Letter is structured as follows. Section [2] intro¬ 
duces the CP-even fraction P+, derives the relations that are used to mea¬ 
sure its value at the -0(3770) resonance, and reviews how knowledge of P+ 
allows non-CP eigenstates to be cleanly employed in the measurement of 
7 with —)■ DK^ decays. Section describes the data set and event 
selection. Sections ii and presents the determination of P+ using CP 


2 


tags, D —it'g lVt+tt” tags and other tags, respectively. In Sect. combi¬ 
nations of the individual sets of results are performed for each signal mode; 
for D —)■ and D —)■ these combinations include the results 

from Ref. |T]. Section gives the conclusions. 


2. Measuring the CP content of a self-conjngate ZJ-meson decay 
and the consequences for the 7 determination with DK^ 

Let the amplitude of a meson decaying to a self-conjugate hnal state 
/ be written as A{D^ —)■ /(x)) = 03 , 6 *®^, where x indicates a particular 
point in the decay phase space and 9x is a CP-conserving strong phase. The 
amplitude is normalised such that 



|^(P°^/(x))|Mx = P(/), 


( 1 ) 


where B{f) is the branching fraction of the decay and V indicates the 
entire phase space. The decay amplitude at x is denoted where x 

indicates the point in phase space reached by applying a CP transformation 
to the hnal-state system at x. CP violation in the charm system is neglected, 
which implies that the decay amplitude at x is equal to the amplitude 
at X. It is useful to dehne the strong phase difference AO^ = 9x — 9x- 

It is possible to express the CP-even fraction in terms of the decay am¬ 
plitudes introduced above. Let the CP eigenstates be \Dcp±) = (|P°) ± 
|P °))/\/2 and consider the decay —)■ / in terms of these states. The total 
CP-even fraction of the inclusive decay is dehned as 


and so 




/xec K/WI^CP+)Pdx 


Lep l(/(x)I^CP+)P + |(/(x)|PcP-)Pdx’ 


Pi = 


+ + 2axax cos A9x dx 


1 + 


Lev 2(«x + ai) dx 


arttx cos A9x dx 


( 2 ) 


( 3 ) 


Note also that the following relation is always true in the absence of CP 
violation: 


ttxttx sin A9x dx = 0 . 


i^x.e'D 


( 4 ) 
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Now consider a quantum-correlated DD system produced in the decay 
of a '0(3770) meson. One of the D mesons in the system decays to / at the 
point X, the other to 5 ^ at y, where in general the phase space of the two 
decays is different. The amplitude of the latter decay is denoted in 

analogy with the terminology used above. 

The amplitude of the '0(3770) —DD fg correlated wavefunction can 
be written [ 6 ] 

•^(/(x)|c/(y)) = ^ . (5) 

The resulting decay probability is then 


^(/(x)|c/(y)) OC 

a^bg + a^by — 2axbyaxby ( cos ^9x cos A 0 j^ -|- sin ^9x sin A 0 y 


( 6 ) 


If both D mesons decay to the same hnal state the probability is divided 
by two to avoid double counting. This formula can be used to determine 
the population of quantum-correlated decays either integrated over all phase 
space or after dividing the phase space into bins. 

The number of ‘double-tagged’ candidates in which one D meson decays 
to / and the other to g, integrating over the phase space of each decay, is 


MU\g) = ZB{f)B{g) 


1 


{2Fl - 1) (2F0 - 1) 


(7) 


where ^ is a normalisation constant common to all decay modes. An impor¬ 
tant special case, considered in Sect. is where the tagging-mode g is a. CP 
eigenstate, and {2FI — 1) reduces to ±1. Section describes an analysis of 
classes of double-tags where this is not the case. 

Alternatively, when the tagging-mode g is a multibody decay, its phase 
space may be divided into bins. Integrating over the phase space of / results 
in the following decay probability in bin i of the phase space of g: 

'P{f\9i)^[ bl + b^- ( 2 FI-I) byby cos A(j)ydy, ( 8 ) 

where Di indicates the phase space encompassed by bin i. In Sect. this 
relation is exploited for the tags D —)■ TTg L^r+Tr”. 

To understand the relevance of the CP-even fraction in the measurement 
of the unitarity-triangle angle 7 consider the decay of a B~ meson to DK~, 
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following which the D meson decays to a self-conjugate hnal state / consisting 
of three or more particles. The amplitude of the B~ decay is a superposition 
of two decay paths: 

A{B-) = A{B- D^K-)A{D^ ^ /) + ^ D^K-)A{W /). (9) 


Following the formalism developed above, the decay amplitude of the 
meson at the point x in the phase space is denoted . The decay 
amplitude of the B~ meson at this point in phase space is therefore 


>f(S-(x)) = A{B- D^K-) 




( 10 ) 


where tb and 6b are respectively the ratio of moduli and the strong phase 
difference between the suppressed and favoured B~ decay amplitudes. The 
resulting decay probability is 


V{B (x)) (X + 2rBaxax cos {6b — 'y + Ox — Gx) 


( 11 ) 


= al + r‘^al + 2rBaxax 


cos{6b — 7 ) cos A 6 * 3 ; — sin(5B — 7 ) sin 


The expression for B~^{x.) is identical except that the sign in front of 7 is re¬ 
versed and X -H- X. The total yield of B^ decays is determined by integrating 
over the entire D phase space: 


Y^ = h^ iP(S^(x))dx 
ixec 

.,2 I fn 17/ 


= 


l + rj^+ {2Fl -l)2rB cos{6b T 7 ) 


( 12 ) 


where is a normalisation constant and Eqs. and have been employed. 
This expression is very similar to that derived in Ref. [7j for the case when 
the D meson decays to a CP eigenstate and is indeed identical in the event 
= 0 or 1. Hence measurements of and observables built from these 
yields pQ, can be used to obtain information on the angle 7 and the other 
parameters of the B^ decay, provided that is known. In Re f. [T] it is 


demonstrated how the effects of mixing, neglected in Eq. lA, may also 

be accommodated. 


3. Data set and event selection 

The data set analysed consists of e+e“ collisions produced by the Cor¬ 
nell Electron Storage Ring (CESR) at a/s = 3.77 GeV corresponding to an 
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Table 1: Z?-meson final states reconstructed in this analysis. 


Type Final states 

Mixed CP 7r~^7i~7r~^7i~, K^K~7r^, Kgj^7r~^7r~ 

CP-even K~^K~, 7r~^7r~, 

CP-odd _ K^7r°, K^u, K^rj, K°r/' 


integrated luminosity of 818 pb“^ and collected with the CLEO-c detector. 
The CLEO-c detector is described in detail elsewhere [S]. Monte Carlo sim¬ 
ulated samples of signal decays are used to estimate selection efficiencies. 
Possible background contributions are determined from a generic sim¬ 

ulated sample corresponding to approximately hfteen times the integrated 
luminosity of the data set. The EVTGEN generator [9] is used to simulate 
the decays. The detector response is modelled using the GEANT software 
package [10]. 

Table [T] lists the P-meson hnal states considered in the analysis. Double¬ 
tag candidates are reconstructed in which one D meson decays into 7 r’'' 7 r“ 7 r’'' 7 r“ 
and the other into a CP eigenstate, or where one D meson decays into 
vr+TT^Tr+TT”, 7 r’'' 7 r“ 7 r° or and the other into one of the mixed-CP 

modes iCgTr+vr” or . The combinations 7 r’'' 7 r“ 7 r’'' 7 r“ vs. 7 r’'' 7 r“ 7 r’'' 7 r“ 

and 7 r’'' 7 r“ 7 r’'' 7 r“ vs. 7 r''' 7 r“ 7 r° are also reconstructed. 

The unstable hnal state particles are reconstructed in the following decay 
modes: 7r° —)■ 77 , ATg — )• tt+tt", uj — ?• tj — )• 77 , 7 — )• Tr+Tr^vr® 

and 7]' —)■ 7 ( 77 ) 7 r+ 7 r“. The 7 r°, Kg, uj, rj and rj' reconstruction procedure is 
identical to that used in Ref. m- 

Final states that do not contain a are fully reconstructed via the beam- 
constrained candidate mass, rubc = \/ s/(dC) — p%/C, where pd is the D- 
candidate momentum, and AE = Ed — \fsl‘2, where Pd is the P-candidate 
energy. The mtc and AE distributions of correctly reconstructed P-meson 
candidates peak at the nominal P° mass and zero, respectively. Neither 
AE nor rubc distributions exhibit any peaking structure for combinatoric 
background. The double-tag yield is determined from counting events in 
signal and sideband regions of rubc after placing requirements on AE mE]- 
[T3] . The selection criteria of candidates involving the modes P —)■ K^K~ and 
P —>■ 7 r''' 7 r“ do not include the cosmic ray muon and radiative Bhabha vetoes 
that are described in Ref. [T] . This is because these sources of background do 
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not contaminate the double-tag sample, and the vetoes are found to perturb 
the selection efficiency of the other D meson in the event. When selecting 
D —>■ Kg7i~^7i~ candidates it is demanded that the Kg decay products form a 
vertex that is signihcantly displaced from the e’''e“ collision point; in contrast, 
for D —)■ 7r’''7r“7r’''7r“ and V 7r~^7r~7r^ candidates the 7r’''7r“ vertex must 
be consistent with originating from the collision point in order to suppress 
contamination from D —>■ Kg7r~^7r~ and D —?■ K^tt^ decays, respectively. 

The double-tag yield determination procedure is identical to that pre¬ 
sented in Refs. mm except for the selections where the signal decay is 
7 r’''7r“7r’''7r“ and the tag decay is K^K~, n~^7r~, or Tr+Tr^Tr+vr", which 

are all dominated by a background from continuum production of light quark- 
antiquark pairs. For these modes an unbinned maximum likelihood £t is per¬ 
formed to the distribution of the average rrihc of the two D-meson candidates. 
The background is modelled with an ARGUS function |T3] and the signal is 
modelled with the sum of two Crystal Ball functions [15] with power-law 
tails on opposite sides. The parameters of the Crystal Ball functions are 
hxed from hts to large samples of simulated data. 

Figures (a) and (b) show the average mbc distributions for CP-tagged 
D —>■ 7r+7r“7r+7r“ candidates, summed over all tag modes that are CP-even 
and CP-odd eigenstates, respectively, where the CP-tag hnal state does not 
contain a meson. Figure]^ shows the average mbc distributions for D 
7r~^7T~7r~^n~, D —)■ 7r’''7r“7r° and D —)■ candidates tagged with D 

Kg7T~^7r~ decays, while Figs. (a)~(c) show the Dalitz-plot distributions of 
the tag decay for these three signal modes. 

Many mesons do not deposit any reconstructible signal in the detector. 
However, double-tag candidates can be fully reconstructed using a missing- 
mass squared (mPgg) technique [IB] for tags containing a single meson. 
Yields are determined from the signal and sideband regions of the mC^g 
distribution. Figure [^(c) shows the mC^g distributions for D —7r+7r“7r+7r“ 
candidates tagged with either a K^tt^ or decay. Figure shows the 
mPgg distributions for D —)■ Tr+Tr^Tr+vr”, D —7r’''7r“7r° and D — 
candidates tagged with D — K^7T^7r~ decays, and Figs. (d)-(f) show the 
corresponding tag-side Dalitz-plot distributions. 

In events where more than one pair of decays is reconstructed an algo¬ 
rithm is applied to select a single double-tag candidate based on the informa¬ 
tion provided by the mbc and AE variables. The particular choice of metric 
varies depending on the category of double tag and is optimised through 
simulation studies. 
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Figure 1: Distributions of D —>• 7r+7r“7r+7r“ candidates tagged by CP-eigenstates. Sub¬ 
figures (a) and (b) show average rrihc distributions for CP-even tags and CP-odd tags 
not involving mesons, respectively. Sub-figure (c) shows the distribution for 

candidates tagged by CP eigenstates that contain a meson. The shaded histogram is 
the estimated peaking background and the vertical dotted lines indicate the signal region. 
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Figure 2: Average Wbc distributions for (a) D —>■ 7r+7r“7r+7r“, (b) D —)• 7r+7r“7r° and (c) 
D —)• iF+iF“7r° candidates tagged by a P —)■ ATgTr+Tr” decay. The shaded histogram is 
the estimated peaking background and the vertical dotted lines indicate the signal region. 


The peaking background estimates are determined from the generic Monte 
Carlo sample of events. For double tags involving a CP mode without 

a meson the peaking backgrounds are found to constitute 5-10% of the 
selected events, and are predominantly from residual D —iCgTr+vr" contam¬ 
ination. The peaking backgrounds for hnal states with a are generally 
larger; for K^tt^ and this contamination amounts to 15-20% of the sig¬ 
nal yield, whereas for K^7r~^7i~ it is ~ 10% of the signal yield. The dominant 
source of peaking background in each case is the equivalent decay containing 
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Figure 3: Dalitz-plot distributions for D —>■ iFg7r+7r“ reconstructed against (a) D —> 
7r+7r“7r+7r“, (b) D —)■ 7r+7r“7r° and (c) D —)■ K^K~'k^, and D —> iF^Tr+Tr” reconstructed 
against (d) D —>■ 7r+7r“7r+7r“, (e) D —)■ 7r+7r“7r° and (f) D —> iF+iF“7r°. The axis labels 
are the invariant-mass squared of the Tr^iFgL pair. 


'400 


o 


100 


,.! ' " 

■ ^ 

(a) : 

- : 

■ t i 

- 

- + 

- 



” + 




- 0.2 0.2 0.6 
<iss [GeVVc^l 




Figure 4: distributions for (a) D —>■ 7r+7r“7r+7r“, (b) D —>■ 7r+7r“7r° and (c) D —)■ 

candidates tagged by a D —>■ iFL7r+7r“ decay. The shaded histogram is the 
estimated peaking background and the vertical dotted lines indicate the signal region. 
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a Kg instead of a K^ meson. The contamination from specific modes in the 
other categories of donble tags is typically 10% or less. The statistical nn- 
certainties on these backgronnd estimates arising from the hnite size of the 
simnlated samples are inclnded in the total statistical uncertainties on the 
signal yields. 

The measured double-tag event yields after background subtraction are 
given in Table 

Table 2: Double-tagged signal yields after background subtraction. Information on the 
entries marked ‘f’? not studied in the current analysis, can be found in Ref. 



+ 

1 

+ 

1 

-1- 

1 

o 

K+K-tt^ 

K+K- 

19.3 ± 6.3 

t 

t 

TT^'K~ 

3.3 ± 8.2 

t 

t 


18.6 ± 5.2 

t 

t 


49.2 ± 10.9 

t 

t 

Klu 

22.0 ± 6.5 

t 

t 


112.8 ± 11.0 

t 

t 


41.0 ± 6.8 

t 

t 

KsVM 

18.8 ± 4.5 

t 

t 

Kgri{7r~^'K~7r^) 

3.1 ± 2.7 

t 

t 

Kli 

9.3 ± 3.3 

t 

t 


217.9 ± 16.8 

289.2 ± 17.6 

52.5 ± 7.8 


485.0 ± 26.3 

291.1 ± 19.2 

78.1 ± 11.1 

+ 

1 

+ 

1 

41.0 ± 16.3 

75.5 ± 15.7 

— 


Knowledge of the single-tag yields of the CP-eigenstate modes is required 
for normalisation purposes. Since the single-tag reconstruction criteria ap¬ 
plied are identical to those employed in Ref. [1], all information on these 
yields is taken from the earlier publication. It is also necessary to know the 
single-tag yield for the decay D —>■ 7 r+ 7 r“ 7 r°. A £t to the mbc distribution 
returns a result of 29998 ± 320 signal candidates, after the subtraction of 
small peaking-background contributions. 

4. Analysis with the CP tags 

The yields of the single and double tags containing a CP eigenstate are 
used as inputs to determine the CP-even fraction, P^^. Following on from 
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Eq. the expected number of observed events, M, where one D meson decays 
to the 7r’''7r“7r’''7r“ hnal state, and the other decays to X, a CP eigenstate 
with eigenvalue r/cp, is given by 


M(47r|X) = 2Njj^B{47i)B{X)e{47i\X) 


1 - r,cp(2Ff - l) 


(13) 


where is the number of DD pairs, i3(47r) and B{X) are the branching 
fractions for the two reconstructed hnal states and £(47r|X) is the efficiency 
of reconstructing such a double tag. The double tag yield is denoted by 
M~ (M+) for CP-even (CP-odd) tags. Experimentally it is advantageous to 
eliminate dependence on the branching fractions and the reconstruction 
efficiency, which can be achieved by normalising by the single-tag yields. The 
yield of single tags, {S~) decaying to a CP-odd (CP-even) eigenstate X, 
is given by 

S(X)^2N^bB(X)s(X), (14) 


where e{X) is the reconstruction efficiency of the single tag. The small 
effects of mixing are eliminated from the measurement by correcting 

the measured single-tag yields 5'^eas ^^ch that = 5'^eas/(1 ~hcp2/p)) where 
Ud = (0.62 ± 0.08)% is the well-known P-mixing parameter [T7]. A further 
correction is applied in the case of the tags K^K~ and tt+tt" because of 
the differing selection requirements for the single and double-tag case, as 
described in Sect. This correction factor is determined by taking the ratio 
of the selection efficiency of the single tag from simulation with the two 
differing selections. It is determined to be 1.15 and 1.10 for the D —)■ K^K~ 
and D —)■ 7r’''7r“ modes, respectively, with an uncertainty of ±0.05. The 
other uncertainties on these single-tag yields are assigned following the same 
procedure described in Ref. pQ. 

For the case of the two CP tags involving a meson a different treatment 
is required, since it is not possible to measure the single-tag yield directly 
for these modes. Following the procedure described in Ref. [T], the effective 
single-tag yield is evaluated using Eq. where the effective single-tag ef- 
hciency e{K^X) is calculated from the ratio of e{ATT\K^X)/e{4:TT), and the 
leading systematic uncertainties are associated with the branching fractions 
and the value used for the effective single-tag efficiency. The effective single¬ 
tag yields are determined to be 21726 ± 3497 and 9124 ± 4105 for and 
AlCJ, respectively. 

Assuming that the reconstruction efficiencies of each D meson are inde¬ 
pendent, then the ratio of the double-tagged and single-tagged yields are 
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Figure 5: 13 —>■ 7r+7r“7r+7r“ results for (a) and (b) N~. In each plot the vertical yellow 
band indicates the value obtained from the combination of all tags. The black portion 
of the uncertainty represents the statistical uncertainty only while the red represents the 
total. 


independent of the branching fraction and reconstruction efficiency of the 
CP tag and This ratio is dehned as N~^ = with an analogous 

expression for N~. The CP-even fraction is then given by 



N+ 

N+ + N-' 


(15) 


The measured values for and N~ for each CP tag are displayed in Fig. 

It can be seen that there is consistency between the individual tags for each 
measurement. The mean value < N~^ >= (5.54 ±0.46) x 10“^ is significantly 
larger than < N~ >= (1.80 ± 0.32) x 10“^, indicating that the 7r+7r“7r+7r“ 
hnal state is predominantly CP even. 

If the acceptance across the phase space of the D —Tr+Tr^Tr+vr" decay 
is not uniform it has the potential to bias the measurement of Using 
simulated data the selection efficiency of individual pions in P —)■ 7r’''7r“7r’''7r“ 
decays is determined in bins of momentum and polar angle with respect to the 
beam direction. The candidates in data are then weighted by the normalised 
efficiency. Each pion is treated independently and an overall weight, typically 
lying within 5-10% of unity, is found by multiplying the individual weights. 
The scaled signal yields are used to re-determine F^^ and the difference 
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between this and the value found without efficiency correction is 0.008, which 
is taken as the systematic uncertainty due to non-uniform acceptance. 

Using the CP tags only, and accounting for the correlations between the 
systematic uncertainties, yields = 0.754 ± 0.031 ± 0.021, where the hrst 
uncertainty is statistical and the second systematic. 


5. Analysis with the D —>■ tags 

For each of the signal samples that are tagged by H —)■ or 

D —)■ AlTt+tt” decays the Dalitz plot of the tag mode is divided into eight 
pairs of symmetric bins by the line , where is the invariant-mass 

squared of the Tr^iFg ^ pair. The bins lying on one side of this line (m^ > m?_) 
are labelled — 1 —)■ —8, and those on the other side 1 —)■ 8. The binning 
dehnition follows the ‘Equal /S.5d BABAR 2008’ scheme of Ref. [18], in which 
the boundaries are chosen according to the strong-phase prediction of a model 
developed by the BaBar collaboration [T9|. The expected distribution of 
entries is symmetric and so the analysis considers the absolute bin number 
|i|, which contains the contents of the pair of bins —i and i. 

Following Eq. the expected population of bin |i| for signal decays with 
Ag7r’''7r“ tags is 


M\i\ = h 


Ki + K_i — (2F+ — l]2ci\/ KiK^i 


(16) 


where h is a normalisation factor specihc to the signal category, is the 
flavour-tagged fraction, being the proportion of Ag7r’''7r“ decays to fall in bin 
i in the case that the mother particle is known to be a meson, and q is 

the cosine of the strong-phase difference between and decays averaged 
in bin i and weighted by the absolute decay rate (a precise dehnition may be 
found in Ref. m- The only difference between the form of this expression 
and the case when the signal decays into a pure CP-even eigenstate [6] is the 
prefactor of (2P+ — 1) in the hnal term. 

Similarly, when the tagging meson decays to K^7i~^7t~ then the number 
of double-tag decays produced in bin |z| is 


M|' I = h' 


K' + KC + 2P, - 1 2c:., K'K'_ 


(17) 


where the primed quantities are now specihc to this case. The reversed sign 
in front of the hnal term rehects the fact that the 77° meson is almost entirely 
a CP-odd eigenstate. 
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The values of Cj and c' within these bins have been measured by the 
CLEO collaboration [TH]. The values of the Ki parameters are taken from 
an analysis of the predictions of various 5-factory models [I9H22] presented 
in Ref. [13], and those of the K[ parameters from measurements performed 
with CLEO-c data [ 23 ] • 

The double-tagged samples are analysed to determine the background- 
subtracted signal yield in each Dalitz-plot bin. The distribution of back¬ 
ground between the different bins is assigned according to its category. Flat 
background is assumed to contribute proportionally to the bin area. Peaking 
backgrounds that occur on the signal side affect the distribution of tag decays 
in lTf+tt phase space according to their nature. For example, in the case 
of 5 —)■ RTsTT® decays that are wrongly reconstructed as 5 —)■ 7r^7r“7r°, the tag 
decay will be in a CP-even state and distributed accordingly. Similarly, the 
distribution of decays that are misreconstructed as 

tags is well understood and modelled appropriately. The distribution of the 
vs. events that contaminate the Tr+Tr^vr+Tr” vs. 

7 r“ selection is determined from data by inverting the iPg veto on the 
signal decay. 

It is also necessary to account for relative bin-to-bin efficiency variations 
in the background-subtracted signal yields. The correction factors are de¬ 
termined from simulation and typically differ ^5% from unity. The signal 
yields in each bin after background subtraction and relative efficiency cor¬ 
rection are shown in Table for tags and in Table for 

tags. 

A log-likelihood £t is performed to the efficiency-corrected signal yields 


residual PoTt+tt 


of each sample, assuming the expected distributions given by Eqs. ^ and [17] 
The £t parameters are the CP-even fraction and the overall normalisation. 
The values of Ki, Ci and c' are also htted, but with their measurement 
uncertainties and correlations imposed with Gaussian constraints. Separate 
hts are performed for the D —iPgTr+Tr” tags, the D —)■ tags, and 

for both samples combined. Fits to large ensembles of simulated experiments 
demonstrate that the returned uncertainties are reliable and that there is no 
signihcant bias in the procedure. All data hts are found to be of good quality. 
The results are plotted in Fig. |^for the D —)■ iCgTr+vr" tags and in Fig. [^for 
the D —)■ K^7r~^7r~ tags. The numerical results for the CP-even fraction are 
given in Table for D —)■ 7r’''7r“7r’''7r“ and in Table for D ^ it 


71 7r° and 


D K+K-ti^ 


The dominant systematic uncertainty is associated with the distribution 
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Table 3: Double-tagged signal yields vs. ifgTr+Tr” after background subtraction in absolute 
bin numbers of the D —>■ /CgTr+Tr” Dalitz plot. The yields are corrected for relative bin- 
to-bin efficiency variations and then scaled so that the totals match the values in Table 



7r'''7r’ 

TT 

'^71 

TT’*' 

71 

'TT^ 

K+K- 

'TT^ 

1 

30.8 

± 

7.0 

29.9 

± 

6.3 

12.6 

± 

4.1 

2 

19.8 

± 

5.3 

19.1 

± 

4.8 

4.6 

± 

2.6 

3 

16.4 

± 

4.5 

27.2 

± 

5.2 

6.9 

± 

2.5 

4 

10.1 

± 

3.4 

18.5 

± 

4.4 

1.6 

± 

1.5 

5 

55.1 

± 

8.1 

96.9 

± 

10.0 

8.4 

± 

3.1 

6 

21.1 

± 

5.1 

31.2 

± 

5.8 

4.4 

± 

2.4 

7 

27.7 

± 

6.0 

34.6 

± 

6.3 

7.7 

± 

2.8 

8 

36.9 

± 

6.8 

31.8 

± 

6.0 

6.2 

± 

2.5 


of the continuum and combinatoric backgrounds. This is assessed by repeat¬ 
ing the measurement with these contributions switched to those found in 
the sidebands of the signal distributions. Shifts of 0.020, 0.019 and 0.037 
are observed for pKK-K° ^ respectively. The uncertainty as¬ 

sociated with the measurement errors on K[, Ct and c' is estimated by 
re-running the £t with these quantities set as fixed parameters and subtract¬ 
ing in quadrature the new £t uncertainty from that obtained with the orig¬ 
inal procedure. This component is found to be 0.013 for D —)■ 7 r+ 7 r“ 7 r+ 7 r“, 
0.010 for D —>■ 7 r+ 7 r“ 7 r° and 0.025 for D —)■ and is accounted 

as a systematic uncertainty in the final results. An uncertainty is evalu¬ 
ated to account for non-uniformities in acceptance across phase space. For 
D —)■ 7 r+ 7 r“ 7 r+ 7 r“ this contribution is calculated with the same procedure as 
in Sect. and found to be 0.002 for the joint Kg i^ 7 r~^ 7 r~ £t. For D —)■ 
and D —)■ the acceptance uncertainties are taken to be 0.001 and 

0.010, respectively, as determined in Ref. pQ. Other sources of bias are eval¬ 
uated to be small. The total systematic uncertainties are given in Tables 
and 1 ^ and are in all cases significantly smaller than the corresponding statis¬ 
tical uncertainties. 

6. Other tags 

The double-tagged yield of 7 r+ 7 r“ 7 r+ 7 r“ vs. 7 r+ 7 r“ 7 r° can be used to de¬ 
termine benehting from the well-measured value of . The ratio 
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Table 4: Double-tagged signal yields vs. after background subtraction in absolute 

bin numbers of the —>■ Dalitz plot. The yields are corrected for relative bin- 

to-bin efficiency variations. 



7r'''7r 

■7r'''7r 

7r~^ 

TT 

71^ 

K+K- 


1 

134.1 

± 

13.9 

89.2 

± 

11.1 

17.3 

± 

6.1 

2 

59.2 

± 

8.9 

32.9 

± 

6.9 

8.8 

± 

4.0 

3 

55.4 

± 

8.7 

31.0 

± 

6.3 

4.1 

± 

3.1 

4 

20.3 

± 

5.8 

7.0 

± 

3.7 

0.1 

± 

1.9 

5 

46.0 

± 

8.7 

6.7 

± 

4.8 

2.1 

± 

3.1 

6 

24.6 

± 

6.2 

14.7 

± 

5.0 

10.0 

± 

3.7 

7 

61.2 

± 

9.0 

46.7 

± 

7.8 

17.6 

± 

4.7 

8 

84.1 

± 

10.8 

62.9 

± 

8.9 

18.1 

± 

5.1 


Table 5: The fit results for the D —>■ if^Tr+Tr tags, where the first uncertainty is 
statistical and the second systematic. The row Xg^Tr+Tr” indicates the configuration 
where the CP-even fraction is a common fit parameter shared between the D —^ Prg7r“'"7r“ 
and D —^ PT^Tr+Tr” samples. 


Tag 

pin 

77g7r’''7r“ 

0.828 ± 0.074 ± 0.014 


0.670 ± 0.057 ± 0.039 


0.737 ± 0.049 ± 0.024 


of double-tag and D —)■ single-tag yields is defined as = 

M( 47 r| 7 r 7 r 7 r°)/S'( 7 r 7 r 7 r°), where a very small correction is applied to the mea¬ 
sured single-tag yield to account for mixing effects. Following from Eq. 
the ratio /N^ removes dependence on the signal branching fraction 

and reconstruction efficiency and is given by 


P^'K'K'K^ 

N+ 



(18) 


which can be rearranged to yield 



yy+ JPTVTTTT^ 

yy'TTTTTrO _ yy'+ _|_ 27 V+ ' 


(19) 
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Absolute bin number Absolute bin number Absolute bin number 


Figure 6: Data (points) and fit results (solid line) in absolute bin numbers for iFg7r+7r“ 
tags vs. (a) D —>■ 7r+7r“7r+7r“, (b) D —i- 7r+7r“7r° and (c) D —i- K^K~'k^. Also shown in 
each case is the expectation if = 0 (dotted line) or = 1 (dashed line). 





Absolute bin number 


Absolute bin number 


Absolute bin number 


Figure 7: Data (points) and fit results (solid line) in absolute bin numbers for 

tags vs. (a) D —>■ 7r+7r“7r+7r“, (b) D —i- 7r+7r“7r° and (c) D —>■ K^K~'k^. Also shown in 

each case is the expectation if 7^+ = 0 (dotted line) or 7^+ = 1 (dashed line). 


The choice of in the denominator of Eq. 1^ is preferred to N~ as it is 
measured with better relative precision. 

Taki^ as input the yields given in Sect. the value of iV+ reported 
in Sect. |^and the hnal result for presented in Sect, [^implies = 

0.695 ± 0.050 ± 0.021, where the uncertainties are statistical and systematic, 
respectively. The main contributions to the systematic uncertainty arise 
from: the measurement of the D —)■ single-tag yield and small violations 

of the efficiency-factorisation ansatz assumed in Eq. f^ the understanding 
of the peaking background component in the sample; and the possible effects 
of non-uniform acceptance. 
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Table 6: The and fit results for the D —>■ itr°7r+7r“ tags, where the first 

uncertainty is statistical and the second systematic. The row lTt+tt” indicates the 
configuration where the CP-even fraction is a common fit parameter shared between the 
D —>■ PlgTr+TT” and D —>■ it'LTr+Tr” samples. 


Tag 


pKKir^ 

Kg7l~^7l~ 

1.03d ± O.OSd ± 0.023 

0.573 ± 0.152 ± O.OdO 

K^TT~^Tr~ 

0.971 ± 0.075 ± 0.033 

0.916 ± 0.181 ± 0.066 


l.Old ± 0.0d5 ± 0.022 

0.73d ± 0.106 ± 0.05d 


The self-tagged yield of vr+Tr^Tr+Tr" vs. 7r’''7r“7r’''7r“ also carries informa¬ 
tion on the value of . This sample is however only used for a consistency 
check, as there are large backgrounds from both the continuum and from 
misidentihcation of D —)■ decays that are a potential source of signif¬ 

icant systematic bias. Furthermore, the predicted yield and measurement 
uncertainty means that the result from analysis of these double tags would 
have low weight in the combined measurement of Using Eq. the 

number of observed self-tagged events is given by 

M (dvr I dvr) = d7^F^’^ (l - , (20) 

where 71 = N^^B{4:7iYe{4:7i\4:7i). The predicted double-tagged yield using 
the value of F^^ obtained from the CP tags is 17±2, which is consistent with 
the measured yield reported in Table 

7. Combination of results 

The results for F^^ from the CP tags, the Kg j^7r~'~7r tags and the 7r'''7r vr® 
tag are summarised in Table They are compatible and are therefore com¬ 
bined, taking account of correlated uncertainties. Correlations arise from 
the non-flat Dalitz plot acceptance between all three measurements and the 
use of N~'~ as an input to both the CP tags and 7r’''7r“7r‘^ tag measurements. 
There is a further small correlation between the results obtained with the 
CP and Tr+TT^TT® tags, associated with the uncertainty on the value of the 
mixing parameter yr). The hnal result is F^ = 0.737 ± 0.028. 

Table [^summarises the results on and obtained with Kg-^7i^7i~ 

tags, together with those determined from CP tags. The iFgL7r’''7r“ mea- 
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Table 7: Results for for each tag category, and combined. When two uncertainties 
are shown, the first is statistical and the second systematic. For the combined result the 
total uncertainty is given. 


Tag 



Fl^ 



CP eigenstates 

0.754 

± 

0.031 

± 

0.021 

^s,L7r+7r- 

0.737 

± 

0.049 

± 

0.024 

7r’''7r“7r‘^ 

0.695 

± 

0.050 

± 

0.021 

Combined 

0 

.737 ± 0 

.028 


surements confirm the results of the earlier analysis. A combination is per¬ 
formed, accounting for the sole source of correlated uncertainties, which is 
that arising from the non-uniform acceptance over the Dalitz plots. Results 
of = 0.973 ± 0.017 and = 0.732 ± 0.055 are obtained. The 

tags improve the relative precision on by 6 % and on 

by 10 %. 


Table 8: Results for F^'^'^° and F^^^° for each tag category, and combined. The CP- 
eigenstate tag results are from Ref. [T]. When two uncertainties are shown, the first is 
statistical and the second systematic. For the combined result the total uncertainty is 
given. 


Tag 

^TTTTTT*^ 


CP eigenstates 

0.968 ± 0.017 ± 0.006 

0.731 ± 0.058 ± 0.021 


1.014 ± 0.045 ± 0.022 

0.734 ± 0.106 ± 0.054 

Combined 

0.973 ± 0.017 

0.732 ± 0.055 


8. Conclusions 

A first measurement has been made of the CP-even fraction of the de¬ 
cay D — 7 r+ 7 r“ 7 r+ 7 r“, exploiting quantum-correlated double-tags involving 
CP-eigenstates, a binned Dalitz-plot analysis of the modes D —)■ ATgLTr+Tr”, 
and D —)■ 7 r+ 7 r“ 7 r° decays. The result, F^ = 0.737 ±0.028, when considered 
alongside the relatively high branching fraction, indicates that this channel is 
a valuable addition to the suite of D decays that can be harnessed for the mea¬ 
surement of the unitarity-triangle angle 7 through the process DK^. 
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The decays D have also been employed as a tag to measure the 

CP contents of the modes D —)■ Tr+Tr^vr® and D —)■ . The results 

conhrm the conclusion of a previous analysis P, based on CP-eigenstate tags, 
and also suggested by earlier amplitude-model studies [231 - I26] . that the CP- 
even content of the 7 r+ 7 r“ 7 r° hnal state is very high, and therefore this decay 
too is a powerful mode for the measurement of 7 . Combining the two sets 
of measurements yields = 0.973 ± 0.017 and = 0.732 ± 0.055. 

Now that their CP-even fractions have been measured, all three decay modes 
may also be used for studies of indirect CP violation and mixing in the 
system [3]. 
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